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INTRODUCTION 


The  adsorption  of  carbon  monoxide  on  nickel  surfaces  has  been  the 
subject  of  a  number  of  experimental  and  theoretical  investigations.  Most 
of  the  experimental  work  and  all  of  the  theoretical  studies  have  involved 
the  low  coverage  regime.  However,  conditions  in  the  high  coverage  regime 
(CO  pressure  >  5  mm)  lead  to  the  reaction  of  CO  with  the  surface  to  form 
nickel  carbonyl:^ 

Ni  +  4  CO  ^  Ni(C0)4(g)  .  (1) 

This  reaction  is  an  intermediate  step  in  the  commercial  purification  of 
2 

nickel.  Furthermore,  it  serves  as  a  prototype  for  most  of  the  technolog¬ 
ically  important  catalytic  surface  reactions,  which  typically  take  place 
under  high  coverage  conditions.  This  system  has  the  advantage  that  dissoci¬ 
ative  chemisorption  and  accumulation  of  product  species  are  not  primary 
sources  of  surface  contamination. 

One  area  of  particular  interest  in  the  nickel  carbonyl ati on 
reaction,  as  well  as  other  catalytic  reactions,  involves  magnetic  effects. 

Enhancement  ("Krinchik  effect")  of  the  reaction  by  external  magnetic  fields 

3  4 

has  been  reported  by  Soviet  workers,  *  but  has  not  been  verified  by 
others. Also,  a  Hedvall  effect  (dependence  of  the  activation  energy  upon 

O 

the  substrate's  magnetic  state)  has  been  observed  for  this  reaction. 

The  detailed  mechanism  of  nickel  carbonyl  formation  has  not  been 

directly  determined.  Although  the  overall  reaction  is  given  by  (1),  it 

almost  certainly  does  not  proceed  directly.  There  is  evidence  that  terminally 

g 

adsorbed  carbonyl  leads  to  NiCO  as  the  first  intermediate  species.  This 
would  presumably  be  followed  by  the  reaction  of  C0(g): 


2 


+  CO  NiCO 

NiCO  +  3  CO  ^  Ni(C0)4  .  (2) 

However,  other  mechanisms  are  possible,  such  as 

N*Surf<“>ads  “  *  Ni(C0)2 

NisurflCOja  *  2  CO  ♦  Ni(CO)^  (3) 

The  fact  that  these  reactions  only  occur  under  high  coverage  con¬ 
ditions  limits  the  sophistication  of  theoretical  treatments.  For  this 
reason  it  is  not  now  feasible  to  address  the  magnetic  state  problems  which 
would  require  a  precise  treatment  of  spin  and  electron  correlation.  We 
instead  have  conducted  the  first  theoretical  examination  of  the  high  coverage 
regime  of  CO  adsorption  on  nickel  in  order  to  gain  insight  into  determining 
the  species  leaving  the  surface  and  details  of  the  reaction  mechanism. 
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II.  DISCUSSION  OF  THE  ELECTRONIC  STRUCTURE  METHOD 

In  problems  such  as  the  one  involving  many  transition-metal  atoms 
and  CO  molecules,  the  application  of  the  usual  first-principles  methods  of 
molecular  quantum  mechanics  becomes  impracticable.  The  large  number  of 
basis  functions,  N,  required  in  this  case  results  in  a  numerical  problem  of 
such  a  magnitude,  O(N^),  that  it  exceeds  present  computer  capacity.  In 
these  cases  one  must  resort  to  more  approximate  methods  which,  while  practical, 
retain  the  prospects  of  yielding  a  realistic  physical  description. 

One  approach  which  has  been  quite  actively  applied  to  electronic 
structure  problems  of  this  magnitude  is  the  approach  of  Anderson. ^ 

This  is  a  semi -empirical  method  which  combines  Hel Imann-Feynman  electrostatic 
forces  and  molecular  one-electron  orbital  energies  to  construct  an  estimate 
of  the  total  electronic  energy  of  a  system.  This  method  has  been  applied 
to  the  study  of  the  structures  of  clusters  of  transition-metal  atoms, to 
the  study  of  chemisorbed  molecules  on  surfaces, and  to  the  study  of 

/14) 

organometall ic  molecules.'  ' 

Procedural ly  the  method  of  Anderson  has  two  steps.  The  first  step 
is  the  calculation  of  pair-wise  repulsive  energies  among  the  atoms 


Er  = 

A>B 


Zr  / 


I  I 

'o/  t/B  ' 


where  is  the  electronic  charge  density  of  atom  A,  Zg  is  the  charge  of  the 
nucleus  of  atom  B,  and  ^  is  the  position  vector  of  atom  A.  The  second  step 
is  the  calculation  of  molecular  one-electron  energy  levels  by  solution,  in 
an  atomic  orbital  basis,  of  the  matrix  wave  equation 
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I  <:j.k  =  0 

J=1 


i  =  1 


l,...,n  ^  N 


In  this  secular  equation,  the  Ej^’s  are  the  molecular  orbital  energies,  the 
S..'s  are  the  overlap  matrix-elements  between  two  atomic  orbitals,  and  the 

^  J 

H..'s  are  the  interaction  matrix-elements  between  two  atomic  orbitals. 

^  *3 

Anderson  defines  the  diagonal  elements  by 

and  the  off-diagonal  (two-center)  elements  by 


■  h  w  'u 

where  is  the  experimental  atomic  ionization  potential  corresponding  to 
removal  from  an  atom  of  an  electron  in  the  i-th  atomic  orbital.  The  factor 
k  has  usually  been  taken  as  2.25  exp(-0.13  D)  where  D  is  the  distance  between 
the  two  atoms. 

The  final  step  in  Anderson's  method  is  to  use  the  molecular  energy 
levels,  E|^,  to  construct  the  energy  term 


n 

^  ^  ^k  ^k 
k=l 

where  pj^  is  the  occupation  number  of  the  k'th  molecular  one-electron  energy 
level.  This  energy  tends  to  be  an  attractive  energy  as  the  atoms  approach 
each  other.  It  is  combined  with  the  previously  defined  repulsive  energy, 
En,  to  obtain  Anderson's  estimate  of  the  total  electronic  energy,  i.e.. 
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The  application  of  this  method  to  large  molecules  has  been  shown 
to  yield  useful  predictions  of  structure,  force  constants,  and  relative 
bond  strengths. ^ There  is  thus  some  empirical  basis  for  its  validity. 

The  evidence  of  its  validity  on  purely  theoretical  grounds  has  been  sparser. 
Anderson  explains  the  physical  reasoning  motivating  his  approach  as  followsi^^^ 
One  first  envisions  constructing  the  molecule  by  superposing  rigid  atoms. 

This  results  in  by  neglecting  in  the  Hel Imann-Feynman  force- formulae  all 
terms  in  the  exact  one-electron  density  not  arising  from  the  superposition  of 
atomic  densities,  and  integrating  the  remaining  force  terms.  The  physical 
effect  which  has  been  neglected  in  the  first  step,  according  to  Anderson,  is 
the  effect  on  the  energy  due  to  the  "tunneling"  of  the  electron  from  one 
atom  to  the  next.  This  effect,  as  is  well  known,  is  expressed  in  the  off- 
diagonal  Hamiltonian  matrix-el ementr,  H..,  connecting  basis  states  describing 

*  J 

electrons  moving  about  alternative  nuclei.  This  is  the  source  for  the 

(15) 

binding  force  between  atoms,'  *  and  is  also  often  referred  to  as  the 
"delocalization"  energy  in  chemistry.  In  the  second  step  of  the  method, 
Anderson  then  identifies  the  delocalization  energy  contribution  to  the  total 
energy  as  the  net  energy  lowering  in  the  orbital  energies  obtained  by  solving 
the  secular  equation.  This  is  subjected  to  the  additional  simplifying 
assumption  that,  for  this  purpose,  field  effects  due  to  neighboring  atoms 
may  be  neglected. 

A  comparison  of  Anderson's  method  with  the  standard  approaches 
is  very  difficult  and  complex.  There  are  a  few  simple  examples,  however,  for 
which  some  contact  can  be  made  with  the  more  conventional  approximations. 
Consider  the  case  of  the  molecular  ion.  This  is  a  one-electron  system. 


0 


Let  us  consider  just  two  basis  states.  Each  basis  state  consists  of  two 
protons  and  one  electron  moving  about  one  or  the  other  proton,  i.e. 


n>  =  N  exp(-|r-R;^l) 

I2>  =  N  expC-ir-Rel 

where  N  is  a  suitable  normalization  factor.  The  secular  equation  is 


In  this  case  ^22  ^  ’  ^12  ^  ^21  ’  ^11  ~  ^22*  lower  energy  is 

^  "  (1+S^2^  ^^ll'^^lZ^ 

The  Hamiltonian  operator  in  this  example  is 


H  =  - : - : - + - : - 

liC-Sfll  lr«el  IWbI 

with  the  property  that 

Hll>=  -I|l> 

H|2>=  -112? 

where  I =-13.6  eV  is  the  first  ionization  potential  of  H  atom.  It  follows 


-I  +<1|  ^ - 11> 


that 
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The  second  term  is  indeed  the  Ej^  term  in  Anderson's  approach.  Therefore 
we  write 

In  reducing  Ej^  one  finds 


where  This  shows  that  ^  0  and  is  thus  indeed  always  a 

repulsive  energy  term.  The  off-diagonal  term  is 


12> 


Since  H^2  ’s  exponentially  dependent  on  the  distance  between  the  atoms,  one 
can  see  that  H.j2  will  give  a  progressively  attractive  energy  lowering  which 
we  identify  as  the  source  of  binding  energy.  Schematically  the  variation  of 
these  energy  terms  with  interatomic  distance,  D,  is 


8 


Anderson's  energy  expression  in  this  case  would  have  been 

^  ^  (1+5^2^ 


where 


=  -I 


^12  "^^12 


One  can  see  that  this  is  almost  identical  to  the  variational  result.  The 
approximation  lies  in  the  approximate  relation 


<1 


-1 


^A“^B I 


2>^  Et,  S 


R  ■"12 


This  shows  that  at  least  in  a  simple  case  one  can  discern,  among  the  energy 
terms  of  a  more  conventional  approximation,  the  physical  effects  which  under¬ 
lie  Anderson's  method,  and  that  they  do  take  on  approximately  the  form  he 
postulates. 
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III.  CO  ADSORPTION  ON  Ni(IOO)  AT  LOW  COVERAGE 

Our  first  studies  examined  the  binding  of  nickel  atoms  in  the  (TOO) 
crystal  face.  A  nine  atom  cluster  was  chosen  for  this  study: 


With  this  structure,  the  central  upper  layer  nickel  atom  modeled  a  typical 

surface  atom  complete  with  its  eight  nearest  neighbors. 

Because  the  orbital  occupation  on  a  nickel  atom  in  the  solid  is 

g 

approximately  d  s  the  repulsions  between  atoms  can  be  modified  by  changing 

(18) 

the  s-orbital  exponent.  Since  one  of  the  energies  of  interest  is  the  energy 
to  remove  a  single  surface  nickel  atom,  the  s-orbital  exponent  was  adjusted 
so  that  the  equilibrium  position  of  the  central  upper  layer  nickel  atom  (as 
it  was  moved  up  or  down  with  all  other  atoms  fixed)  lay  in  the  plane  with 
the  other  four  upper  layer  atoms.  The  same  exponent  was  used  for  all  nine 
atoms.  Figure  1  shows  this  procedure.  With  an  exponent  of  -  2.2,  the 
repulsions  are  too  small  and  the  central  atom  contracts  below  the  surface. 
Conversely,  at  1.9,  the  atom  pops  out  of  the  surface  plane.  A  value  of 
1.97  was  found  to  be  nearly  optimal.^  The  binding  energy  for  removing  the 
central  nickel  atom  (Nig  -►  Nig  +  Ni)was  found  to  be  0.5  eV.  This  does  not 

agree  well  with  the  experimental  bulk  value  of  5.6  eV.  This  problem  is 

'  ■  ■ 

^  An  error  in  tabulation  led  us  to  believe  that  this  value  was  closer  to  the 
optimum  than  2.0.  Either  of  these  values  is  reasonable,  but  the  final 
result  may  be  slightly  ^  2.  This  point  will  be  corrected  in  future  studies. 
It  should  be  noted  that  the  difference  in  binding  energies  calculated  with 
1.97  and  2.0  is  only  .03  eV.  See  also  the  results  in  Section  VI. 
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partly  linked  to  problems  in  parameterization  end  in  occupation  numbers. 

In  generating  the  potential  curve  the  d-orbital  energy  was  kept  fixed  at 
the  bulk  value  -10.00  eV.  Since  the  s-orbital  energy  was  -7.63  eV,  this 
resulted  in  a  d^^  nickel  atom  (with  total  energy  -100  eV)  being  removed. 

If  the  gas  phase  nickel  parameters  were  used  for  the  single  atom,  its  total 
energy  would  be  only  -77  eV,  a  difference  of  23  eV!  Thus,  even  switching 
parameters  smoothly  along  the  curve  would  not  lead  to  a  reasonable  result. 
This  problem  is  less  noticeable  in  the  removal  of  bonded  species  where 
orbital  mixing  occurs  or  the  d^^  configuration  becomes  plausible.  For 
example,  the  tetrahedral  geometry  of  Ni(CO)^  can  be  interpreted  as  nickel 
in  a  d  configuration  with  its  vacant  s  and  p  orbitals  sp  (tetrahedral ly) 
hybridized  to  participate  in  bonding  to  the  CO  lone  pairs.  However,  the 
binding  energy  is  known,  both  theoretically^^  and  experimentally^^,  to  be  a 
function  of  the  cluster  size.  Thus,  it  is  not  surprising  that  the  nine  atom 
calculation  does  not  yield  the  bulk  value. 

Standard  parameters  were  used  for  C  and  0.  The  set  of  parameters 
is  given  in  Table  1.  The  bond  length  for  CO  was  never  varied.  A  value  of 

O 

1.15  A  was  used.  For  the  nickel  surfaces  a  nearest  neighbor  distance  of 

o  o 

2.4918  A  (lattice  edge  =  3.52388  A)  was  employed  throughout. 

Using  these  parameters  and  bond  lengths,  the  adsorption  of  CO  at 
various  sites  on  the  (100)  surface  was  examined,  (See  the  left-hand  part  of 
Figure  2.*)  The  optimal  height  above  the  surface  was  determined  for  each 
site.  The  curves  in  the  figure  are  plotted  as  a  function  of  the  resultant 

Figure  2  will  be  discussed  further  in  Section  V. 
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Table  1.  Orbital 

Exponents  and  Energies 

Orbital 

Exp 

e(eV) 

Ni 

4s 

1,97 

-7.63 

4p 

1,5 

-4.00 

3d* 

5.75 

-10.00 

11 

2.0 

C 

2s 

1.61 

-16.6 

2p 

1.57 

-11.3, 

-11.2^ 

0 

2s 

2.25 

-28.5 

2p 

2.23 

-13.6 

Contraction  coefficients  were  0,5683  and  0,62922 


^  Inadvertently,  a  value  of  -11,2  eV  was  used  for  the  low  coverage 
(100)  studies.  This  oversight  was  corrected  in  subsequent  cal¬ 
culations. 
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Figure  2.  CO  Adsorption  Equilibria  at  Low  Coverage. 
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Ni-C  bond  length(s).  Even  though  the  CO  comes  down  to  only  1.1  A  above  the 
nickel  surface  for  the  4-fold  site,  it  still  has  the  longest  Ni-C  bond 

O 

lengths,  2.06  A.  The  4-fold  site  for  the  chosen  9  atom  cluster  was  modeled 
by  CO  approaching  from  the  underside.  The  bridge  site  used  here  was  one  of 
the  ones  on  the  top  side.  These  optimizations  were  necessary  because  keeping 
the  CO  at  a  constant  distance  from  the  surface  favors  one  site  over  the 
others. 

The  net  energy  'aH)  required  for  the  removal  of  a  surface  atom 
(from  the  same  nickel  cluster)  by  a  terminally  adsorbed  carbonyl  molecule 
(NigCO  Nig  +  NiC0)was  calculated  to  be  1.1  eV. 

Using  the  same  parameters,  we  conducted  a  series  of  calculations 
of  the  total  energies  (stabilities)  of  Ni(C0)p  for  n  =  1-4.  All  spin  con¬ 
figurations  used  were  closed  shell.  The  Ni-C  bond  length  used  was  that  of 

O 

nickel  carbonyl,  1.87  A.  The  geometries  for  n  =  2  and  3  were  those  of  a 

(20) 

similar  study  using  the  method  .  Our  results  are  shown  in  Table  2. 

These  values  are  required  for  the  energetic  analyses  in  Section  IV-D.  They 
show  Ni(CO)^  to  be  the  most  stable,  with  stability  decreasing  as  the  number 
of  bonded  CO's  decreases  (in  contrast  to  the  calculations  ).  This  seems 
only  natural  since  the  final  product  of  reaction  (1)  is  the  tetracarbonyl . 
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Table  2,  Binding  Energies  of  Nickel  Carbonyl  Species 


Species 

Total  Energy  (eV) 

Binding  (eV)"^ 

Binding/co 

AB 

Ni  +  CO 

-278.182 

— 

— 

— 

Ni(CO) 

-279.996 

1.814 

1.814  eV 

1.814 

Ni(C0)2 

-459.651 

3.287 

1.644 

1.473 

Ni(C0)3 

-639.509 

4.963 

1.654 

1.676 

Ni(C0)4 

-819.105 

6.377 

1.594 

1.414 

^  relative  to  Ni  +  nCO 
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IV. _ CO  ADSORPTION  ON  Ni(IOO)  AT  HIGH  COVERAGE 

A,  Geometry  Model 

,  v(21) 

LEED  studies  of  carbon  monoxide  adsorption  on  nickel  (100)  have 

determined  that  for  coverages  (e)  up  to  50%  the  adsorbing  molecules  form  a 

simple  ^2x2  (Figure  3a).  By  e  =  .61  the  adsorbed  layer  has 

shifted  to  a  hexagonal  pattern  which  is  compressed  by  3%  (Figure  3c).  The 

maximum  compression  observed  corresponds  to  a  coverage  of  0.69  monolayers 

15  2 

(which  is  equivalent  to  1.10  x  10  molecules/cm  ).  For  intermediate  coverages 

(50-60%)  some  disorder  was  indicated. 

With  these  facts  in  mind»  a  model  for  high  coverage  was  devised. 

Figure  3  illustrates  the  relation  between  the  symmetries  of  the  substrate  and 

overlayer.  The  LEED  studies  do  not  determine  whether  the  €2^2  overlayer 

structure  (Figure  3a)  should  be  registered  over  terminal,  bridge  (shown),  or 

four-fold  bonding  sites.  The  intersections  of  the  lines  indicate  CO  molecules, 

and  nickel  atoms  are  indicated  ::v  dots.  Our  studies  on  NigCO  indicate  that 

location  at  the  four-fold  sites  is  preferable.  The  elongated  hexagonal 

pattern  in  Figure  3b  is  obtained  by  shifting  the  CO  positions  in  every  other 
+  ° 

column  of  3a  by  -  1.762  A.  If  this  pattern  is  then  compressed  horizontally 
(see  arrows)  by  one  nickel  column  in  every  six,  the  60%  coverage  pattern  in 
(c)  results.  This  structure  is  compressed  by  about  5%  from  a  regular  hexagon. 
When  the  compression  is  one  of  every  four  columns  of  nickel  atoms,  the 
resultant  coverage  is  67%.  These  two  coverage  levels  are  close  enough  to  the 
critical  coverages  found  experimentally  to  suggest  that  the  periodic  regis¬ 
tration  which  defines  the  models  (despite  their  general  misregistered  character) 
is  significant,  and  that  the  models  are  reasonable. 
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B.  Results 

The  most  sophisticated  clusters  for  the  high  (60%)  coverage  cal¬ 
culations  consisted  of  a  single  layer  of  nine  nickel  atoms  overlaid  by 
three  CO  molecules  in  a  nearly  equilateral  arrangement  (see  Figure  4).  Due 
to  periodicity,  it  was  not  necessary  to  vary  the  x-coordinates  more  than 

O 

1.175  A  in  the  course  of  seeking  the  optimum  position  of  the  CO's  relative 
to  the  substrate.  This  was  an  important  factor  in  minimizing  edge  effects. 

In  most  cases,  energies  from  two  calculations  (labeled  by  J^eft  and  nght) 
were  required  to  obtain  an  energy  representative  of  the  hexagonal  pattern, 
since  the  environments  of  adsorbed  molecules  in  columns  A  and  B  of  Figure  3c 
are  different  and  there  are  twice  as  many  B  columns  as  A  columns.  Arrange¬ 
ments  were  labeled  by  the  site  of  the  central  hexagonal  molecule,  as  indicated 
in  Figure  3.  Table  3  gives  the  adsorption  energies  obtained  for  the  ^erminal- 
_bridge-,  and  four- fold  centered  arrangements  and  the  three  intermediate 
arrangements.  In  these  calculations,  the  height  of  the  adsorbed  molecules 
above  the  nickel  surface  was  variable.  Since  it  was  not  feasible  to  actually 
optimize  each  z-coordinate,  each  was  approximately  optimally  predetermined  by 
a  model  based  on  the  low  coverage  optimized  distances.  The  contours  shown 
in  Figure  5  indicate  the  values  used.  The  arrangement  of  concentric  square 
contours  corresponds  to  a  3-D  description  in  terms  of  square  pyramidal  peaks 
and  holes  at  the  terminal  and  four-fold  sites,  respectively.  This  flexibility 
in  the  cluster  model  was  necessary  because  using  a  fixed  distance  would 
unfairly  favor  certain  sites  over  others.  The  results  of  Table  3  show  that 
the  bridge-centered  arrangement  was  the  most  stable,  but  that  due  to  the 
forced  misregistration,  all  of  the  energies  were  fairly  similar. 
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Figure  4.  x-  and  y-Coordinates  of  High  (60%) 
Coverage  Model  Calculations. 


20 


Table  3. 

Adsorption  Energies  of  60%  Coverage 

Model  Calculations 

Location 

Total  Energy  (eV) 

Average  Adsorption 
Energy  (eV)/C0 

B 

-1439.783 

2.82 

B/4 

-1439.099 

2.59 

4 

-1438.896 

2.53 

4/T 

-1439.474 

2.72 

T 

-1438.959 

2.55 

T/B 

-1439.181 

2.62 

The  question  remains  as  to  how  the  coverage  switches  from  4-fold 
(centered)  ^2x2  to  bridge-centered  hexagonal.  Suppose  that  Figure  3a  had 
been  drawn  4-fold  centered.  Then  the  shifting  process  (to  obtain  a 
hexagonal  pattern)  shown  in  (b)  would  have  resulted  in  half  of  the  adsorbed 
molecules  being  in  terminal  sites  (which  are  the  least  favorable  at  low 
coverage).  Thus  rearrangement  in  this  fashion  is  unlikely  for  energetic 
reasons.  Of  course,  without  compression  a  shift  to  the  bridge-centered 
hexagonal  would  not  occur  either.  The  compression  occurs  because  the 
stability  gained  by  the  increased  number  of  absorptive  bonds  more  than  makes 
up  for  the  decrease  in  binding  per  bond  when  a  less  optimal  site  is  used. 
After  the  4-fold  site,  the  bridge  site  is  the  next  best,  so  the  result  is  a 
combination  of  4-fold  and  bridge  sites.  The  hexagonal  pattern  minimizes  the 
intermolecular  repulsions  between  adsorbed  molecules. 

C.  Theoretical  Coverage  Limit 

In  principle,  it  should  be  possible  to  estimate  the  absolute 
coverage  ratio  theoretically.  Some  complications  do  arise  from  the  use  of 
a  finite  cluster  model,  however.  Consider  that  for  an  infinite  number  of 
CO  molecules  in  a  two-dimensional  hexagonal  arrangement,  the  number  of 
nearest  neighbor  interactions  per  molecule  is  three.  However,  our  three- 
molecule  model  hcs  only  one  per  molecule.  A  seven-molecule  cluster  would 
still  have  less  than  two.  Thus,  the  proper  amount  of  repulsion  for  the 
infinite  (macroscopic  or  experimental)  case  is  not  closely  modeled  by  small 
clusters.  This  would  lead  to  an  overestimate  of  stability. 

When  ^he  nickel  substrate  is  truncated, the  cluster  is  left  with 
"dangling  bonds",  i.e.,  the  nickels  are  coordinated  by  fewer  nearest 
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neighbors  than  usual.  This  applies  to  truncation  in  all  three  dimensions 
(ix,  iy,  and  -z,  excepting,  of  course,  the  intended  surface).  This  might 
also  enhance  the  coverage  calculation  by  overestimating  the  binding  energy 
per  CO  molecule.  Thus,  unless  adjustments  are  made,  the  theoretical  coverage 
value  will  probably  be  an  upper  limit.  However,  since  the  interpretation  of 
LEED  data  is  often  subject  to  discrepancies  of  a  factor  of  two,  a  theoretical 
bound,  even  an  approximate  one,  would  be  of  interest. 

Since  the  favored  geometry  at  60%  coverage  was  the  "bridge  centered" 
one  we  start  from  that  point.  The  study  consisted  of  (arbitrarily)  further 
compressing  the  lateral  dimension  (separation)  of  the  CO's,  starting  from  a 
sheet  of  nine  nickel  atoms  and  three  CO  molecules.  The  energies  for 
coverages  of  50%,  60%,  67%,  and  100%  were  computed.  The  67%  geometry  was 
included  because  it  was  a  symmetric  point  and  it  was  near  the  experimental 
absolute  coverage  value.  The  stability  is  calculated  by  determining  the 
binding  energy  per  CO  molecule  and  multiplying  it  by  the  coverage  (shown  in 
column  3  of  Table  4).  Its  units  are  therefore  eV  per  nickel  atom,  and  the 
predicted  absolute  coverage  corresponds  to  the  maximum  in  the  stability. 

These  calculations  found  the  100%  case  to  be  the  most  stable!  Further  exam¬ 
ination  revealed  that  the  intermolecular  repulsions  between  neighboring  CO's 
were  entirely  too  small,  so  that  even  trying  to  compensate  for  the  missing 
nearest  neighbor  repulsions  would  be  fruitless. 

Therefore  a  correction  procedure  was  devised.  The  original 
binding  energies  are  given  in  the  first  column  of  Table  4.  For  the  distorted 
hexagonal  structure  at  the  60%  coverage  geometry,  there  are  two  nearest 
neighbors  (nn)  plus  one  almost  as  near.  At  100%  the  structure  will  be  ^2x2 
with  two  nn  and  two  next  nn  (nnn)  interactions  per  CO.  The  (nnn)  analogs 
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Table  4.  Calculated  and  Adjusted  Stabilities  as  a  Function  of  Coverage 


Coverage 

Binding  Energy/ 

CO 

Stability 

eV/Ni 

Adjusted  Stability 
eV/Ni 

B-hex 

50% 

2.988 

1.494 

1.480 

B-hex 

60% 

2.820 

1.692 

1.649 

B-hex 

67% 

2.562 

1.708 

1.609 

B-C2X2 


2.963 


2.963 


0.557 
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in  the  hexagonal  structure  are  the  almost  nn  and  a  single  nnn.  At  50% 
hexagonal  coverage,  the  pair  of  nn  distances  are  actually  larger  than  the 

O 

single  ("almost")  nn  and  the  nnn  distance  is  >  7  A.  The  repulsive  part  of 
the  interaction  potential  between  CO  molecules  (fit  to  experiment)^ 

^  6.7829  X  10^  eV 
“^repulsive  ~  ^12 

can  be  used  to  correct  the  stabilities,  (The  entire  potential,  including  the 
attractive  term. 


^  -68.6  eV 

’attractive  3 
r 

could  be  used.)  The  correction  per  CO  molecule  according  to  the  appropriate 
neighbors  described  above  has  been  calculated  and  weighted  by  the  coverage  to 
give  the  adjusted  stability  values  shown  in  Table  4.  These  crude  corrections 
indicate  the  absolute  coverage  level  for  the  bridged  models  examined  to  be 
closer  to  60%  than  100%, 

In  the  future,  ab-initio  calculations  of  the  CO-CO  interaction 
should  be  carried  out  to  be  sure  that  the  potential  above  is  not  the  result 
of  an  erroneous  assumption  of  the  coverage  level.  Subsequently,  the  CO 
parameters  (orbital  exponents)  used  in  the  semi -empirical  calculations  should 
be  modified  to  reproduce  this  interaction  more  reliably.  This  reparameter¬ 
ization  would  be  expected  to  decrease  the  adsorption  energies,  bringing  them 
closer  to  agreement  with  experiment.  It  is  possible  that  the  ordering  of 
site  preference  would  also  change. 
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D.  Energy  Analyses 

It  is  possible  to  analyze  certain  aspects  of  the  formation  of  nickel 
carbonyl  from  an  energetics  viewpoint.  Consider  a  nickel  surface  with  a  given 
amount  of  CO  already  adsorbed  upon  it  and  an  influx  of  additional  CO 

What  products  are  likely?  Since  we  have  calculated  the  average  adsorption 
energy  (per  CO)  at  several  coverages  and  total  energies  of  the  various  nickel 
carbonyl  species  we  can  examine  the  relative  energetics  of  several  processes. 
Suppose  we  use  the  adjusted  adsorption  energies  (as  a  function  of  coverage) 
given  in  Table  5.  We  assume  that  n  and  m  are  large  enough  that  removal  of 
one  nickel  or  addition  of  carbonyl  does  not  affect  the  coverage  significantly. 
Figure  6  shows  the  relative  energetics  of  some  of  the  competing  processes,  as 
a  function  of  coverage.  To  make  the  curves  basically  horizontal  we  do  not 
include  the  energy  due  to  addition  adsorption  of  CO  as  m(0)  increases,  since 
it  involves  a  large  quantity  which  all  of  the  processes  have  in  common.  In 
Figure  6  we  see  that  at  low  coverage,  the  most  favorable  process  is  further 
adsorption.  Moreover,  in  contrast  to  nickel  carbonyl  formation,  we  expect 
no  activation  energy  for  this  process.  However,  as  the  absolute  coverage 
limit  is  approached,  the  energetic  advantage  of  this  product  channel 
decreases.  Eventually  the  limit  is  reached  so  that  for  every  CO  that  adsorbs 
another  desorbs.  At  this  point,  the  most  favorable  product  is  Ni(C0)4. 

Unfortunately,  this  graph  does  not  tell  us  which  nickel  species 
first  desorbs.  We  cannot  tell  from  it  the  magnitude  of  the  probable  barriers 
to  reaction.  From  our  calculations  we  can  conclude  that  even  at  higher  cover¬ 
ages  there  is  probably  not  a  significant  barrier  to  additional  adsorption 
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Table 

5.  Adjusted  Adsorption  Energies 
as  a  Function  of  Coverage  for 

Bridge-Centered  Model  Sites 

0 

AE(eV) 

0.0 

-3.12 

0,5 

-2.96 

0.6 

-2.75 

0.67 

-2.41 

0,7 

-2.0 

0.75 

-1.5 

0.8 

-0.5 

>0.8 

-0.0 

/ 


Nin_i+3(C0)  +  Ni(CO) 
Nijj_j^+2(CO)  +  Ni(CO) 

Nin_i+(CO)  +  Ni(CO), 
^in-1  Ni(CO)^ 


0.5 


1.0 


Coverage 


Figure  6,  Relative  Energetics  of  Competing  Processes  as 
a  Function  of  Coverage. 


29 


because  of  the  flatness  of  the  potential  energy  surface  due  to  misregistra¬ 
tion.  We  may  also  examine  the  optimal  geometry  for  clues. 

Let  us  assume  that  at  67%  coverage  that  the  bridge  centered 
arrangement  is  still  the  most  favored.  Consider  the  situation  of  an  addi¬ 
tional  CO  molecule  approaching  the  surface.  Actually,  it  first  nears  the 
layer  of  already-adsorbed  carbonyls.  The  most  favorable  landing  sites  on  the 
nickel  substrate  will  presumably  be  determined  by  the  holes  in  the  overlayer, 
rather  than  by  direct  interaction  of  the  incoming  carbonyl  with  the  nickel 
surface.  We  thus  superimpose  a  "landing  pattern"  of  probable  interaction 
sites  on  Figure  3d,  shown  by  dotted  lines  in  Figure  7.  We  now  indicate  the 
adsorbed  carbonyls  by  open  circles. 

This  pattern  is  another  hexagonal  pattern,  but  it  is  rotated  -v  30° 
from  the  carbonyl  lattice.  We  note  that  for  most  of  the  more  probable 
landing  positions  (particularly  at  the  hexagonal  corners),  the  incoming 
carbonyl  will  land  near  a  nickel  atom  which  is  interacting  with  a  nearby 
(already-adsorbed)  carbonyl.  This  suggests  the  possibility  that  Ni(C0)2  might 
be  the  candidate  desorbing  species.  These  molecules  have  been  indicated  by 
drawing  in  stick  bonds.  There  are  only  two  unique  corner  sites.  One  has  one 
nickel  dicarbonyl  possibility  and  the  other  has  two.  In  these  high  coverage 
calculations  we  have  not  allowed  the  surface  nickels  to  move.  It  is  possible 
that  the  adsorption  of  carbonyl  would  lengthen  the  surface-nickel-atom-to- 
bulk  bond  lengths,  which  would  presumably  enhance  the  possibility  of  reactions. 

Examination  of  high  coverage  on  other  surfaces  is  of  interest 
because  of  recent  experimental  findings  which  reveal  that  due  to  some  unknown 
selectivity,  etching  of  (100)  and  (110)  surfaces  leads  exclusively  to  (111) 
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facets,  thus,  the  steady  state  nickel  carbonyl  formation  rates  observed 
starting  with  various  surface  orientations  have  actually  been  only  the  rate 
on  (111).  Further  experimental  investigations  to  measure  the  rates  on  non- 
(111)  faces  will  be  difficult  because  they  will  presumably  be  concerned  with 
the  early  stages  of  reaction,  when  the  well-known  problems  of  surface 
cleanliness  and  preparation  are  particularly  troublesome.  On  the  other 
hand,  the  difficulty  of  theoretical  studies  is  not  affected  by  the  choice 
of  crystal  orientation.  Thus,  theory  may  be  able  to  aid  in  the  determination 
of  the  relative  ordering  of  rates  for  the  various  surfaces  and,  thereby,  the 
cause  of  facetting.  It  is  possible  that  examination  of  "landing  sites"  will 
lead  to  such  explanations. 
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V; _ PRELIMINARY  CALCULATIONS  ON  THE  (110)  AND  Mil)  FACES  OF  NICKEL 

In  light  of  the  preceding  discussion,  some  preliminary  calculations 
of  CO  adsorption  on  the  (110)  and  (111)  nickel  surfaces  were  performed.  The 
previous  nickel  s-orbital  exponent  [optimized  for  the  (100)  surface]  was 
used  for  preliminary  calculations  on  the  (111)  and  (110)  surfaces.  Surface 
atoms  in  these  orientations  have  different  numbers  of  neighbors  than  the  (100) 
face.  The  (111)  surface  atoms  have  nine  rather  than  eight  nearest  neighbors, 
six  in  the  surface  plane  and  three  in  the  next  layer. ^  The  (110)  surface  is 
furrowed  in  comparison  to  the  (100)  and  (111)  faces.  It  is  composed  of 
alternating  upper  and  lower  rows  of  atoms,  such  that  those  in  the  upper  row 
have  7  nearest  neighbors,  while  those  in  the  lower  (almost  internal)  row  have 
11. 

Binding  (adsorption)  of  CO  at  various  sites  on  small  nickel  clusters 
modeling  the  (100),  (110)  and  (111)  faces  is  illustrated  (as  a  function  of 
Ni-CO  bond  length  in  Figure  2.  The  curves  have  been  shifted  so  that  CO 
adsorbed  at  each  bridge  site  is  the  zero  of  energy.  CO  on  Ni(lll)  has  basi¬ 
cally  three  types  of  bonding  sites:  terminal  (T),  bridged  (B),  and  three¬ 
fold  (3)  sites  with  either  a  hole  or  an  atom  underneath  (in  the  second  layer). 
The  Ni-C  bond  lengths  were  optimized  at  these  sites  (only  the  3-fold-atom 
site  was  used  since  the  3-fold-hole  site  was  assumed  to  be  similar  because 

o 

second  layer  location  is  'v-  3.4  A  away).  As  for  the  (100)  face,  the  highest 
coordination  (3-fold)  site  was  found  to  be  preferred.  Again  the  terminal 


^  The  additional  repulsion  which  results  indicates  that  (ideally)  the  exponent 
might  need  to  be  increased  for  nickel  atoms  on  this  surface.  However,  a 
small  difference  here  may  not  be  significant  at  this  point. 
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site  was  the  least  favored,  and  the  Ni-C  bond  length  was  similar  to  that  for 
the  Ni(lOO)  terminal  site. 

On  Ni(llO)  there  are  narrow  (B)  and  wide  (W)  bridge  sites,  a 
terminal  site  (T),  and  a  rectangular  (R)  site.  Here  also,  terminal  bonding 

O 

near  1.8  A  was  less  ideal  than  the  other  sites.  The  most  favorable  site  was 
the  wide  bridge  site.  The  rectangular  site  was  intermediate  (between  the 
two  bridge  sites)  energetically. 

These  sites  should  be  examined  more  closely.  The  rectangular  site 
is  actually  a  5-fold  site  since  it  has  terminal  bonding  to  the  Ni  atom  in  the 

O 

trough  row.  The  bond  distances  to  the  four  corner  nickels  are  2.27  A,  and 

O 

the  terminal  bond  is  1.95  A.  The  wide  bridge  site  has  Ni-C  bond  lengths  of 

O 

1.96  A  to  the  peak  row  nickel  atoms  and  2.44  to  the  trough  row  nickel  atoms, 
and  could  thus  be  considered  to  have  4-fold  coordination.  The  C-(peak/ 
trough)  Ni  bonds  for  the  short  bridge  site  vary  more  (1.91/3.22)  and  do  not 
indicate  4-fold  coordination.  For  comparison  the  second  layer  nickel  atom 
for  the  (100)  4-fold  site  are  2.86  A  from  the  CO,  and  the  two  for  the  (100) 

O 

bridge  site  are  >  3.4  A  away, mostly  because  the  distance  between  layers  is 
0 

0.5  A  larger  for  (100),  The  five-fold  site  is  probably  not  the  most  favorable 
site  because  it  has  mainly  terminal  bond  character  (less  favorable)  and  the 

bond  distances  to  the  rectangular  corner  nickels  are  stretched  in  comparison  to 
the  (100)  4-fold  case.  Moreover,  the  wide  bridge  site  bears  some  similarity 
to  the  (favored)  4-fold  site  on  (100):  the  4-fold  site  can  be  considered  as 
a  "double  wide-bridge  site"  because  of  its  diagonal  Ni-Ni  distance.  Subsequent 
calculations  will  use  readjusted  nickel  and  carbonyl  parameters  to  determine 
these  single  carbonyl  adsorption  geometries,  and  the  character  of  the  bonding 
at  each  site  will  be  analyzed  in  detail. 
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VK _ EMPIRICAL  MODEL 

This  last  section  describes  an  empirical  modeling  program  which  has 
been  assembled  for  the  purpose  of  locating  reaction  pathways  and  transition 
states  for  the  various  processes  relevant  to  nickel  carbonyl  formation.  The 
economical  advantages  of  this  model  include  being  able  to  treat  much  larger 
clusters  and  to  calculate  many  more  points  on  the  potential  energy  surface. 
Within  this  model,  CO  is  treated  as  an  atom.  Its  parameters  are  adjusted  to 
reproduce  the  adsorption  energy,  frequency,  and  equilibrium  geometry  of  CO 
on  nickel . 

In  our  first  application  of  this  model,  a  5  layer  cluster  of  46 
atoms  represented  the  nickel  surface.  When  the  top  layer  was  allowed  to 
relax,  the  surface  equilibrium  position  (inter-layer  distance)  was  found  to 
be  3%  larger  than  the  bulk  separation.  This  result  implies  that  the  use  of 
nickel  exponents  (in  the  semi-empirical  calculations)  which  correspond  to 
longer  surface  nickel  bonds  lengths  may  actually  be  more  realistic  than  that 
of  those  which  reproduce  the  bulk  characteristics. 

A  second  preliminary  calculation  modeled  the  desorption  of  NiCO 
from  the  nickel  (100)  surface.  The  system  was  treated  in  a  three-body  LEPS 
(London-Eyring-Polanyi-Sato)  formulation,  in  which  the  components  were 
a  single  Ni  atom,  the  (n-l)-atom  nickel  cluster,  and  the  CO  molecule.  The 
interaction  parameters  (in  this  case)  were  fit  to  experimental  values  for 
nickel  and  adsorbed  CO,  and  ab  initio  values  for  NiCO.  The  model  predicted  an 
intrinsic  barrier  of  0.7  kcal/mol  resulting  in  an  activation  energy  (relative 
to  adsorbed  CO)  of  127.6  kcal/mole.  This  quantity  is  the  total  of  the 
binding  energy  of  a  nickel  atom  to  the  surface  plus  the  adsorption  energy  of 
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CO  on  the  nickel  surface  minus  the  molecular  binding  energy  of  NiCO  [relative 
to  Ni(g)  +  C0(g)]  plus  the  intrinsic  barrier  height.  These  results  correspond, 
of  course,  to  the  low  coverage  case.  As  a  first  guess,  one  would  expect  that 
the  activation  energy  (relative  to  free  CO)  under  high  coverage  conditions 
would  be  primarily  unchanged,  and  that  the  barrier  relative  to  adsorbed  would 
approach  this  value  as  the  binding  per  additional  CO  tended  to  zero.  Thus, 
these  preliminary  calculations  indicate  a  rather  substantial  activation  energy 
requirement  for  nickel  carbonyl  formation  via  the  monocarbonyl  intermediate. 

This  model  will  be  used  in  future  studies  to  further  examine  the 
dicarbonyl  and  other  arrangements  suggested  by  landing  site  arguments  and  for 
subsequent  transition  state  investigations. 
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